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ABSTRACT 


An error probability analysis is done for a DFT based, M-ary frequency-shift 
keying (MFSK) communications system employing fast frequency-hopped spread 
spectrum signals. A linear combination procedure referred to as noise-normalization is 
employed at the receiver to minimize the effects of partial-band interference, which is 
modeled as additive Gaussian noise. The performance of the receiver is studied as a 
function of signal Doppler shift and type of windowing function used in the DFT. 

The use of fast frequency-hopped spread spectrum is found to improve the 
performance of the DFT based receiver in all but the most severe cases of Doppler shift. 
The use of a non-uniform window (i.e.. a Hamming window) to improve receiver 
performance 1s effective only in the presence of large Doppler shifts. The amount of 
Doppler shift necessary to warrant the use of a non-uniform window depends on the 
amount of jamming noise power at the receiver, but is relatively insensitive to the 
frequency-hop rate used. In general, in the absence of any information concerning the 
mature of the received signal, a non-uniform window should be used because the 
performance degradation experienced at small Doppler shifts is insignificant compared 


to the performance enhancement gained at large Doppler shifts. 
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I. INTRODUCTION 


A. NON-UNIFORM WINDOWING OF THE DISCRETE FOURIER TRANSFORM 

The use of the discrete Fourier transform (DFT) to 
noncoherently detect signals employing M-ary frequency-shift 
keying (MFSK) modulation has gained popularity recently due to 
advances made in high speed digital systems. At the current 
level of technology, however, the utility of the DFT based 
receiver remains limited due to the speed of real-time DFT 
algorithms. Even the performance enhancement obtained by using 
the fast Fourier transform (FFT) is not sufficient to overcome 
this limitation. However, certain applications - namely fast 
frequency-hopped spread spectrum systems - can effectively 
utilize the DFT receiver because lower data rates frequently 
correspond to higher immunity to narrowband interference. 
eet. Isp. 20147) 

Problems arise in using the DFT based receiver if the 
transmitted frequency does not match the assigned frequency of 
the DFT bin. This situation results when ee ghee. frequency 
is different than that required by the detector, as when the 
Signal experiences a Doppler shift prior to detection. This 
leads to a situation where transmitted signal energy falsely 
contributes to an incorrect DFT frequency sample: a phenomena 


known as self-induced crosstalk [Ref. l:p. 2015]. This problem 


is not specific to the DFT based receiver. What is unique 
about this type of receiver is its ability to minimize self- 
induced crosstalk through the use of non-uniform windows. An 
example of a DFT based detector experiencing self-induced 
crosstalk is presented in Fig. 1.1. Figure 1.1(a) depicts the 
Situation where the received frequency is perfectly matched to 
the first output bin of the detector, while Fig. 1.1(b) shows 
the results of slightly modifying the received Signal's 
frequency. As can be seen, the magnitude of the bin 
corresponding to the signal is reduced, while the magnitudes 
of neighboring bins which represent other frequencies are 
Significantly larger. Large amounts of frequency deviation can 
easily lead to unacceptable bit error rates. 

One method that is used to minimize the effects of self- 
induced crosstalk involves the use of nonuniform windowing; 
that is, time-sampled data is not equally weighted. The 
frequency response of typical nonuniform windows have broader 
main lobes than do rectangular windows (uniform weighting) but 
exhibit sharper rolloff characteristics. The sharp rollotit 
characteristics of nonuniform windows partially compensate for 
the effects of self-induced crosstalk, but in some cases the 
existence of a broader main lobe creates distortion where it 
would not otherwise exist [Ref. l:p. 2015]. Nonetheless, 
nonuniform windowing can be used successfully to combat large 
doppler shifts in low data rate environments as demonstrated 


in | BigenawrZar 
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Figure 1.1: Outputs of a DFT (a) without a doppler shifted 
Signal, and (b) with a doppler shifted signal. The adverse 
effects of a doppler shifted signal can severely degrade 
communications. 


B. FACTORS AFFECTING RECEIVER PERFORMANCE 
1. Fast Frequency-Hopping and Partial-Band Interference 
When the carrier of the transmitted frequency is 
changed in an apparently random manner during the transmission 
of a signal, the overall effect is to spread the energy of the 
Signal over a large bandwidth. This spreading technique is 
called frequency-hopping and is illustrated in Fig. 1.2. Two 
types of frequency-hopping are considered: fast and slow 


frequency-hopping. Slow frequency-hopping systems change, or 


hop, the carrier once per data symbol, while fast frequency- 
hopping systems hop the carrier more than once per data 
symbol. The higher the hopping rate, the more the frequency 
Spectrum approximates that of true spread spectrum systems 


utilizing pseudonoise (PN) codes [Ref. 3: pp. 10-14]. 


Carrier Frequency 








Time 


Figure 1.2: An example of hopping the carrier frequency of 
a signal. 


The intent of a hostile jammer is to disrupt 
communications. In order to effectively jam a frequency-hopped 
Signal, a smart jammer may employ partial-band interference; 
although, this type of interference can also be caused by 
other unintended narrowband sources. In utilizing partial-band 


interference, the jammer concentrates all of the jamming 


Piergy into a fraction of the total bandwidth and randomly 
hops this signal through the spectrum in the same fashion as 
the signal employing frequency-hopping. By doing this, the 
frequency hopped signal cannot always avoid the portion of the 
spectrum that is jammed; the jammer must be dealt with via the 
fmececiver. 

The assumptions involved in the analysis to follow are 
that the jammer is limited in power; that the interference can 
be modeled as additive Gaussian noise; and that the 
interference, when present, is present in all detection 
channels with probability y. Hence, the probability that 
interference is not present in any branch is l-y. Thus, y 
represents the fraction of the spread bandwidth that is 
jammed. If the partial-band interference is assumed to have a 
power spectral density of N,/2 over the entire spread 
bandwidth, the amount of interference present at the receiver 
due to partial-band interference is y'N,/2 with probability y 
and 0 with probability l-y. [{Ref. 4:p. 3] 

In addition to partial-band interference, the system is 
assumed to be corrupted by additive white Gaussian noise with 
a power spectral density of N,/2. Accordingly, the power 
Spectral density of the total noise at the receiver input is 
No/2+y'N;/2 with probability y and N,/2 with probability 1-y. 
If the receiver equivalent noise bandwidth is W Hz, then the 
total noise power at the receiver is NW with probability 1l-y 


and is (y7'N,+Ny)W with probability y. [Ref. 4:p. 3] 


2. Fading 

In addition to partial-band interference, the 
Signal is assumed to be further degraded by fading. Fading is 
caused by components of the signal arriving over different 
paths which combine to form the total detected signal. Due to 
the different path lengths traversed by the various signal 
components, the combined signal amplitude may be greater than 
or less than the expected amplitude without fading due to 
constructive and destructive interference. Fading severely 
degrades the performance of a real world communications system 
and must be taken into account in any plausible receiver 

design and analysis. 
For this analysis, the signal is assumed to experience 
Slow, frequency non-selective, Rician fading. Slow fading 
refers to the fact that the amplitude of the signal is assumed 
to remain constant at least over the duration of one hop, 
while frequency non-selective implies that all frequency 
components of the signal experience identical fading. Rician 
fading is characterized as having a diffuse signal component 
and a direct signal component. Assumme that the amplitude of 
the received signal in a Rician fading channel is av2. Then 
the probabilty density function of the random variable A is 


given by 








fylayetie PF 1,(-2Ey ula), (1) 


iijement(s) is the umit stepefunction,»~I" is the power»in the 
direct component of the signal, 20, is the power in the 
Guanusescomponent, and i,(*) is the modified Bessel function 
of order zero. As can be seen, (1) reduces to the well known 


Rayleigh distribution when the value of the direct component 


momzero. [Ref. 5] 


C. THE NOISE-NORMALIZED DETECTOR 
A noise-normalized detector is a linear receiver that uses 
the received noise power to normalize the output of each 
branch before a decision is made regarding which signal is 
present. This type of receiver is combined with fast frequency 
hopping to improve overall system immunity to partial-band 
interference. The noise-normalization is implemented before 
the hop receptions are combined to form the detection 
Statistics. This type of detection scheme tends to lessen the 
influence on the overall decision statistic of each hop when 
interference is present, while increasing the influence of 
each hop when no interference exists [Ref. 4:p. 2]. Because 
the output statistics of each branch of the receiver are 
independent of one another, closed form solutions for the 
receiver performance can be obtained. 
This thesis studies the effects of nonuniform 
windowing on a noise-normalized DFT based receiver where the 
Signal experiences partial-band interference, Doppler shift, 


and Rician fading. The effects of fast-frequency hopping the 


carrier frequency are andlyzeq-and the resules for Varioucsaa. 
rates discussed. As in [1] and [2], a statistical analysis is 
done to determine the amount of frequency deviation necessary 
to warrant the use of a nonuniform window over a rectangular 


window. 


II. ANALYSIS OF THE NOISE-NORMALIZED DFT BASED DETECTOR 


The receiver under consideration is shown in Fig. 2.1. The 
incoming signal hop is first brought back to baseband and then 
sampled at the Nyquist rate. The DFT of the signal is computed 
using the appropriate windowing function and the magnitude of 
each output bin is Squared. Hence, in the absence of both a 
Doppler shift and a nonuniform window, this receiver is 
functionally equivalent to a conventional MFSK receiver 
mipovying quadratic detectors. At the same time, the total 
noise power at the receiver is measured, modified 
appropriately by the windowing function of choice, and 
inverted to form the noise-normalization term. Each bin output 
is multiplied by the noise-normalization term. The decision 
statistics are then obtained by summing over all hops 
GCOmprising the transmitted symbol. Finally, the detection 
Statistics are compared and the largest one used to determine 
the estimate of the symbol sent. Errors result when the 
detection statistic in a channel with no signal is greater 
than the detection statistic of the signal channel. A detailed 


analysis of each part of this design follows. 
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Figure 2.1: The DFT based noise-normalized receiver. 
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A. DETERMINATION OF THE PROBABILITY DENSITY FUNCTION FOR EACH 
DFT OUTPUT BIN 
Pome we bow eNestecelver in Fig. 2.1 every fT, 


seconds is of the form 


r(t) =af/2cos[2n(f,+£,+f') t]+n’(t), (2) 
where O<t<T,, @ iS sa Rician random variable having the 
probability density function of (1), avV2 is the Rician 
eer buced signal amplitude, f. is the carrier frequency 
ieeped every T, seconds, fas the intended symbol frequency, 
f' is a Doppler shift, and n'‘(t) is zero mean, additive 
Gaussian noise. The signal frequencies are chosen to be 
SatEnogonal during the time interval O<t<T, and are obtained 


using the relation 


I 
gill lea (3) 
L'p Dy, 


where m'=1,2,..,M and a toed neocon mr epresentingm™rne 
Spacing between each of the M signals. To simplify the signal 


representation, the following substitutions are made: 


qé1+A-(m’-1) (4) 


and 


eaf’T, (5) 


so that after the low-pass filter (2) becomes 


joe 


r'(t) =ay2cos(2={2*©) ©) sn(e), (6) 
h 


where € is the fraction of Doppler shift with respect to the 
bandwidth of the baseband information signal, g is an integer 
related to the signal frequency £., and™a( 2) isco ef ae ma 
noise. The discrete form of the received signal is obtained by 
using a sampling rate of N/T,, where N is the number of 


samples taken during the symbol reception: 
r (p) =aV2cos [244LSP] +n(p), (7) 


where p=0,1,..,N-1 and is the discrete signal sample. [Ref. 
Z:pp. 7-95 
The output of each bin is a complex quantity and can be 


separated into real and imaginary parts: 
Yn (K|a) =Yom(K/@) +7 ¥ tm (kla) » (8) 


where k, which ranges from 0 to N-1l, is an integer 
representing the DFT output bin and Y (kla), Y,,(k|a), and 
Y ,(k|a) are respectively the DFT output, the real pare of@eme 
DFT output, and the imaginary part of the DFT output for 
channel k and frequency hop m conditioned on the Rician 
distributed random variable a. The real and imaginary parts of 


the DFT ane calculated wsepamarey 


1,3 


Y,,(k|a) > Dee r(p)w(p)cos [ SRP | 


r-1 2™ (Qt€) D 2m Kp 
-y , [V2acos [——4 —“* ] +n(p) ]w(p) cos[—= =] , 


(9) 


and 


Vz, (kla) =}, x (p) w(p) sin [S222] (10) 


="? [y2acos [SEZ P j sn(p) ] w(p) sin{ 2], 





where w(p) represents the windowing function. Assuming slow 
Rician fading so that the amplitude of the signal remains 
constant during each hop interval, the conditional probability 
eemoity function for the random variable X,, is 
— (Xm + Wien) 
Xm 20.4 ua ( kml? km ) (11) 
0 


2 S 2 ‘ 
are 0, 





fy 


km 





ta (Xiem|@) = 


where X,, is the magnitude of bin k and hop m given by 


one] 12 (UR Ey) 


=/Yp,(k\a)*+Y,,(kla)?, 


Peete Che Magnitude Of the expected value of bin «k and 


(12) 


frequency hop m given by 


Him= |Xiem| (13) 


= Vem (K]4) “+¥rm (Ja, 
and o,* is the variance of bin k given by 
0,2==E(|¥,(kla) ~¥, (kay [7] . (14) 
ger. ip. 2016] 


3 


The mean of xX, for a fixed amplitude a a5 ecaciiy 
determined by observing that the only random quantity is Wp) 
which is zero-mean. This leads to the conclusion that the mean 
1s simply the output of the DFT without noise. The real and 


imaginary parts of the mean conditioned on a are 





- : 2mp(q+é) 2mpk 
Yam (Ka) =av2)0 w(p) cos (~*~ — ) cos ( ay ) (15) 





4M, , 
and 
Y,, (kla) =ay2y w(p)cos Aha) sin ( LIES 
m Dp=0 N N . (16) 
SM rm 


An expression for Bh conditioned on a is then determined by 
substituting (15) and (16) into (13) and squaring both sides: 
=m +3 
=a? (m'2+m'4) (17) 
ee ae 
where $8, is the magnitude squared of the mean of the DFT 
output without fading when the signal power is normalized to 


Unity ,ehat cy 


_ -1 2m p (qt+e) 2npk 
Bim=2 ( Dene Wi D) COS ame rant Ne 


(18) 
til ipa cos (SEPLG*E) sin ( SERS) ]2) 


[Rer.w22pe., 10S) 


The variance of X,, is obtained by substituting (9), (10), 


(15), and (16) into (14) and sineliiyeene 
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02-2 EL | ¥qq(Kla) ~Poq KAT +5 (Yam (Kl) ~YoqTRTAY) 1 
e Zin KE 
= ae w(p) n(p) cos (=—*) (19) 


ITN w(p) n(p) sin( FE) |?) 
The magnitude squared of (19) becomes 


= El (Yo w(p)n(p) cos ( SEP) )? 
bie wip) n(p) sin ( 2EXB) | 
=E| yas Dale w (Dp) w(s) cos (=2°P ) cos ( 27 XS) n(p) n(s) 
als Dries 
yee es hAGe) Ql Sp) Selig a ) sin (=—— a lniGaacs) |, 


(20) 





where s is an integer having the same range of values as p. 
Because the expected value operator is a linear operator and 
mies Only random quantities are the noise functions, (20) 


Simplifies to 


2s = yw w(s) cos (3B) 


ein 


> yw s) sin ( a 


BEES} tapi nta. 


OS 





(21) 
Say) ( 
The autocorrelation function is defined as 


Rin(P, S) =E[n(p)n(s)], (22) 

which, because the noise functions represent white noise, is 
Rp (PD, S) =WN,5 (p-s) , (23) 

where 6(:) is the Dirac delta function, W represents the 


equivalent noise bandwidth of the low-pass filter at the front 


kd 


of the receiver, and N. 1s the power spectral densityeereen. 
total noise seen at the receiver [Ref. 4]. The quantity N, 
equals either N, when jamming is not present or NotN,/y when 
jamming is present. The expected value needed to solve (21) is 


obtained by substituting (23) into (22): 


E([n(p)n(s)]=WN,6(p-s) . (24) 


The final expression for the variance Of Xis obtamees 


by substituting (24) into (21) and simplifying: 





_ WN, = =F 2nmkKp 2nukKs 
0,7=— = y..., ¥(p) w(s) cos COs One) 








WN 2 Z : ‘ 
+BY Fg wows) sin AIP) ain ZEEE OC) ay 
mis ye 2, 2UKP, ..:-2, 2nkp | 
= 5 ~ W (p) [cos ‘ree asm ale aia ] 





_ WN, Be 
serge, 7 U2) « 


The variance in each bin is affected by both the power 
spectral density of the external noise and the choice of 
windowing function used in the DFT. 

The probability density function for the output of each 
bin is found by substituting (17) into (11) and integrating 


over the range of the Rician distributed amplitude: 
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cee =| Fy ja (Xpmla) £, (a) da 
~ (Xpn+a2B pp) 
Se 
a ae 0 2 
Ox 0. 
-(a2+T) r 
a 20. a 
=é Hy (eels (26) 
Se Of 
2 
1, Xm, I 
‘im Ss ( gee G 2 











By making the substitutions I,(x)=J)(jx), and by defining 


Pin aot 

20,° 260, 
2 2 

= B im? ¢ +O, 


2q 2 
20,°0 5 


07 4 





(27) 


and using the integral (from [7]) 








oo 2.2 - ( 2 
[ Bae Boe Wetec) oy (2) Cxc— fie eee (28) 
: 2G 20? 


the probability density function is 


f Spe op hod A i 
km 20,70 7 
Z 2 2 
FG, OR 
2 


40? ae . o 4 Xm y B ie 
I, (—=——>) 


sen = 
(29) 








The next step involves substituting (27) into (29) which gives 


1 


2 2 
=Ke.d 7 -) -o . XiemP xm? ¢° +120 ,* 


fy (rex 0 : 20,70 7 20,76 5" (B pnp? +97) 
BimO ¢ +9 (30) 
2 2 
Bm ¢ +9 x 


Equation (30) simplifies to 
2 
y ( Xjen* Bian 
Xen Bn? ro? 7 (_Xeol'VBan y (32) 


—_——_—— CO 
B im ra “ake B em £° ion 


= 
2 


IG) = 


In order to clarify the notation, the following substitutions 





are made: 
op bho -+0,°, (32) 
and 
a, al) ae (33) 
With these substitutions, (31) 18S seen to be Rician 
distributed 
— (Xian? +O tan”) 
_ Xm 20 kn? X km km 34 
EE Se) rt woe Owe ( 
km km 


and reduces to the result obtained in [1}) for the case of a 


nonfading channel (20,4=0 and T=1). [Ref. 2:pp. 11-12] 
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B. ANALYSIS OF THE NOISE-NORMALIZED RECEIVER 
1. Determination of the Probability Density Function for 
the Detection Statistic 
Poemadeed aimee s il, the random variable Zz, is 


Given by 





Z em= j (35) 


where Z,, is the random variable for bin k and frequency hop 
m. Because all the random variables X, are independent, the 
transformation of random variables is accomplished by using 


the relation (from [8]) 


fy (Xkmioy ) 


36 
9" (Xpmioy ) _— 


fy \ Zum) = 


where 


g'(x,,) 1s the derivative of g(x,,) with respext to x,, and 
[ee che real root of (35). The denominator in (36) is 


found by substituting (35) into (37) and differentiating: 


We, 


2Xim (0) 


/ = 
g (Ximiay ) 7 - 2 


(38) 





The transformation of random variables for the noise- 
normalized detection statistic is determined by substituting 


(38) INtOmesiboe 


iN ales 
Lo q §Zkm) = Zi Thee V Fx Zn) (2 


The probability density function of 2. is obtained by 


Substituting (34) into (39): 


3 a (0,2 Z pnt @ em) 


Fp (2j,) = Cay Ten (etn — 
2 km © Km 202 : 2 
O km O km 





The actual detection statistic needed for the 


comparison is Z,, which is defined as 
Z.=y) Pye, (41) 


where L is the number of frequency hops per symbol. Since some 
hops are jammed while others are not, the probability density 
function of Z, conditioned on the number of jammed hops is 
Simply the multiple convolution of the probability density 
functions of the jammed and unjammed hops. For instance, if a 


symbol has three hops and two of them are jammed, the 
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Peobability density function of the detection statistic is 


determined as follows: 
Ey (Zy) = Lz, ( Zeq) FY) Fz, (Zug) APY) Fy, (Zp) amet, = (42) 


where ® represents convolution. In general, the probability 
density function of Z, conditioned on 1 out of L hops jammed 


is 
ee ( ae = eal) jammed) & “@) ep. (a unjammed) & (L-1) ; (43) 


Emere 2c represents ac -fold convolution. [Ref. 4:p. 7] 
2. Determining the Probability of Bit Error 
Once the random variable transformation is 
accomplished, the next step involves calculating the 
probability of bit error. An error will occur if the power in 
a channel other than the signal channel is greater than the 
power in the signal channel. Mathematically, this is expressed 


as 
oo 2; 
Pr (Z,<Z,) ={ £,,(z5) f ie, Ue) Claekn (44) 


where i,j=l,...,M and i#j7; and the signal is assumed to be in 
branch 2.[{Ref. 4: p. 7). The double integral in (44) 
determines the probability that the power present in the 
Signal channel is less than the power present in a noise 


channel. 
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Because the distribution of power in each noise 
channel is not identical when Doppler shift is present, the 
comparison in (44) must be carried out for each noise channel 
and the results combined to give the total probability of 
error. An upper bound to the total conditional probability of 


error is found using the union bound: 


where Ps., is the probability of symbol error given I hops are 
jammed with the signal present in channel i. [Ref. l:p. 2017] 

Because the distribution of power in each bin varies 
with the location of the signal, the amount of Doppler, and 
the relative direction of the frequency shift, an average must 
be taken over all signal channels and direction of Doppler 


shift: 


Ps, (€) ey, ei Ps,,(+€) TES eo), (46) 


where Ps,(€) is the total probability of symbol error given J 
hops are jammed conditioned on the signal experiencing a 
Doppler shift of e« and Ps,,(-) is the probability of symbol 
error given 1 hops are jammed with the signal present in 
channel i conditioned on the signal experiencing a Doppler 
Shift of (+é€) or (-€). The union bound on the conditional 
probability of symbol error is calculated by combining (45) 


and (46): 


Z2 


A 
IPS ((C.) errr) Die ist [Pr(Z,<Z,|+e) +Pr(Z;<Z;|-€)]. (47) 


The total conditional probability of symbol error assuming l 


hops are jammed is calculated by combining (44) with (47): 


4 oo 24 


eo Zz 
+f Fosj-e(23|-€) f £7, \-2 \2q| ©) dz ,dz,] : 


feet. lip. 2017) 
The total probability of symbol error for the noise- 
normalizing receiver in the presence of partial-band 


interference is 


Ps(e)=)>, ( ma le) 7 PSE) (49) 


where Ps(é) 1s the total probability of symbol error [Ref. 
4:p. 6]. To determine the worst case probability of symbol 
error, (49) is evaluated for different values of y, the 
results compared, and the highest conditional probability of 
symbol error recorded. The worst case probability of bit error 


Ss : 


_M/2 
Pb(€) =" Ps(€), (50) 


where Pb(€) is the probability of bit error. 
An exact solution for the conditional probability of 
symbol error is calculated by evaluating the probability of 


not making an error (i.e., the probability that the power in 
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all noise channels is less than the power in a signal 


channel): 


Pr (Zee) =f fz, (23) fo 'f,, (23) dz,dz,. (51) 


The total conditional probability of not making an error is 
found by the following relation: 

PS;, (€) =] [j-1 Pr (Z;<Z;|e) 

Fei 7 (52) 
z yen f £2, (2ile) f f(z leazaen, 

where Ps. ,(€) is the probability of not making an error given 
that the signal in channel i with 1 hops jammed and the signal 
experiencing a Doppler shift of é€. The probability of making 
an error given that the signal is in channel i with J hops 


jammed is then given by 


Ps,,(€) =1-Ps;,(€) - (53) 
which is then used in (46), (47), (48), (49), and (50) to 
obtain the exact probability of bit error. Due to the enormous 
amount of computing power needed to numerically integrate (51) 
and the desire to compare the results with those obtained in 
[1] and [2], the conditional probability of error based on the 


union bound was chosen for this thesis. 
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III. NUMERICAL PROCEDURE 


A. PARAMETERIZATION OF VARIABLES 

In order to compare the performance of the noise- 
normalizing receiver with a conventional DFT/MFSK receiver, 
the number of points used in the DFT, the cut-off frequency of 
the low pass filter, and the spacing between the M signals are 
chosen to be consistent with those used in [1] and [2]. The 


following definitions result from adopting these conventions: 


N44M, (54) 

WA2M, (55) 
and 

A 42, (56) 


where N is the number of signal samples per hop, M is the 
medulation Order, W is the effective noise bandwidth, and 4, 
is an integer representing the message frequency spacing. 
penopetame, 7), Ofmone seconderesudts from sampling the 
mecervea Signal at the Nyquist rate. Whenethis value of T, and 
ey aLle Substituted into (3), it is determined that f, 
consists only of odd values ranging from 1 to 2M-1. 


Mieencrmnore, Ehe Spacangect Emequencies in the DFT is 
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4M (57) 


where Af is the frequency resolution of the DFT [Ref. 9]. 
Because of these choices one empty bin exists between each 
possible DFT signal output. 

A parameter of concern is the energy-per-bit to thermal 
noise power spectral density ratio defined as E,/N,, where E, 
is the transmitted energy per bit. The signal-to-noise ratio, 
which is a measure of the ratio of signal power to noise 


power, is defined as 
+ 
Ni ————__— (58) 


For an M-ary system, symbols are sent instead of bits. The 
power per symbol is the sum of the power in the direct 


component and the power in the diffuse component: 


FE=(20 -4]-) T,, (59) 
where E. is the transmitted energy per symbol and T, 1s the 


symbol duration time related to the hop duration time by 
FIO IT. (60) 


The energy-per-symbol is related to the energy-per-bit by 
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Eee log, (™) . (61) 
All the variables used in this thesis are expressed as per-bit 
quantities; therefore, a relationship between the energy-per- 
hop and the energy-per-bit is needed. Because each symbol is 


comprised of L bits, 





E,=—, (62) 


Muere EF is the energy-per-hop. A relationship between the 
energy-per-bit to noise power spectral density ratio and the 


Signal-to-noise ratio is obtained by combining (58) with (59) 


3 ae (63) 
N, log, (™) 


Seeper choice OL the symbol! duration time (1i.e., T,=1) and the 
equivalent noise bandwidth (1.e., W=log,(M)) gives E,/N,=SNR. 
In general, however, these two ratios are not equal and must 
be kept separate in the receiver analysis. 

An expression relating the energy-per-bit to noise power 
spectral density ratio to the signal-to-noise ratio is 


obtained by combining (58), (60), and (63) and simplifying: 


Ep _ (W)(L)(T,) T+20,° (64) 


‘Ny, log, (™) NW 


For the case of a receiver operating in a jammed environment, 
N, is replaced with NotN,/y when the hop is jammed and (64) 


becomes 
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5, 1 220) | ei aaa, 
Ny 1+N;/(¥Ny)° W(N,+N;/y) log,(M 


(65) 
The relationship between the energy-per-bit to noise power 
spectral density and the signal-to-noise ratio is a function 
of the effective noise bandwidth of the system, the number of 
frequency hops per symbol, the hop duration time, and the 
level of M-ary signalling used. 

Another parameter for analysis is the direct to diffuse 


ratio 





(66) 


The following expressions for the direct and the diffuse 
Signal power components result from combining (65) and (66) 


and simplifying: 


E 
Oo 2 67 
i. 2. (eae) T’, (on 
and 
pe= 2 (_DID (68) 


Jin elie 


In doing this analysis, the values of E,/Np, N;/No, y, W, T, 
and DTD are assumed and used to determine the values of IF* and 


2 
OO. e 
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B. CONVOLUTING THE PROBABILITY DENSITY FUNCTIONS 

In the case of fast frequency-hopping, the probability 
density function of the total detection statistic is the L- 
fold convolution of the probability density functions of the 
detection statistic of each hop. As shown in (49), the total 
probability of error represents all possible combinations of 
jammed and unjammed hops. Based on this, a total of 2M(L+1) L- 
fold convolutions are required for each possible signal branch 
and for each combination of jamming power and fraction of 
bandwidth jammed that are studied. These convolutions are not 
required in the case of slow hopping because the total 
detection statistic 1s equal to the detection statistic of the 
Single hop. These multiple convolutions are generated by 
taking advantage of the properties of the Laplace transform. 


The Laplace transform of the probability density function 


LE Lz im (he 4 =F y(n) (s) 


oo 7 (n) (69) 
“|, E yim Can te peas a 


where n=1,2 and n=1 represents the condition of a jammed hop 
and n=2 represents the condition of an unjammed hop. The 
Laplace transform of interest is obtained by substituting (40) 


into (69) 
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es (0,2 (n) +0 97!) 


Zim 
P10 = rs Boy ¢ im OV Zi i Pe 
im 0 20 2 (n) Gpcne (70) 
5 -@; 2 (n) _ 2) ee pe = ~ 
a 0, 20; 2(n) He 20 oz X sm Ow Zim ) dzin 
20,,2(”) 0 0 a(n) am 


Equation (70) 1S evaluated by defining 


ig. (og 
ps2 im Ox : (71) 
go. 2{n) 
im 
and 
3.2 
Co Sie a (72) 
20g. 


2 


and making the substitutions u°=z, and I)(x)=J)(jx), and using 


the integral (from [{10]}) 





Bes 
fe C™tug, (Bu) du-—+_e 4@, (73) 
0 ee 
The Laplace transform of interest is 
P ; =n '™ g 
We cme (74) 
Py tn) (s) An) = ’ 
where 
A =2 07 ‘™ s+o,? (75) 


[Ref. 4:p. 6]. The probability density Lunction™ ot 92 aa. 
determined by using the fact that convolution in the time 


domain translates to multiplication in the s domain: 


5)10) 


Pye, (sin, (s)*. (76) 


The probability density function of Z, is given by 


£2 (2;) =27 0 (Fm (8) )7 (Fz (s)) 27], (77) 


where L is the number of hops per symbol, 1 is the number of 
jammed hops, and gY' is the inverse Laplace transform 
operation. In general, (77) cannot be evaluated in closed form 
and numerical subroutines are used to determine the inverse 
Laplace transform at each point. 

For the case of all hops jammed or unjammed, (77) reduces 


to 


(n) 
yim (eae 


“wee” 


=2" [ (Fz (Ss) ) 7] 


oo n (78) 
=i (Fi (s))te8" ds. 
0 im 


The function in the s domain in (78) is obtained by raising 


(74) to the L power 


Pe ee ee) 20 im s+o,? (79) 
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Equation (79) is factored as 
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je cp Lb. 205,20") 
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The needed inverse Laplace transform is obtained by combining 


the transform property that 


Ot F(st+a) ] =e? F(z) , (81) 


the inverse Laplace transform relation (from [11]) 


b Loa 
#3 ((4) +e #] =(2) 2 T. (2ybz), (82) 
with the substitutions 
Ge 
303m’ (83) 


and 


pe LM im Ox (84) 


After simplifying, the resulting closed form expression for 


this special case is 


2 Se 0 cc “9x +O (A) 85 
(2740;"7) 2 


(n) ~ re), 0 
Q 20 sn x , [zQ ny) | 


2 (n) 2(n 
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where I, is the modified) Bessel funeriom of (ender ert yand 
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Ss a (86) 


Meer. 4:p. 7}. 


C. NUMERICAL INTEGRATION TECHNIQUES 

A double numerical integration is needed to evaluate (48), 
which can be very time consuming. By taking advantage of the 
properties of the Laplace transform, the double numerical 
integration is reduced to a single numerical integration. The 


amount of computing time needed is greatly reduced since 
%3 _co-l 1 
[7 f2,(ap dzi-8" (Fz, (8) (87) 


and (76) can be divided by s as needed to carry out the inner 
numerical integration. 

The use of the integration technique illustrated in (87) 
is the primary reason a union bound solution is calculated 
instead of an exact solution. As seen in (44), the union bound 
solution requires implementation of (87) only once. The exact 
solution, as seen in (51), requires using (87) M-1 times 
Since, for the DFT based receiver, the bins that ideally 
contain no signal will in general have different probability 
density functions due to Doppler shift and non-uniform 
windowing. Although using (87) to integrate the probability 
density function significantly reduces computing time compared 


to direct numerical integration, the amount of computing 
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resources needed to calculate the exact solution is much more 
than that needed to determine the union bound solution. The 
fact that only a limited amount of these resources was 
available to generate results justified the decision to 
evaluate the union bound solution. 

The remaining integration is accomplished using Simpson's 
rule, which is chosen because of its accuracy as compared with 
the trapezoidal rule. Figure 3.1 shows how Simpson's rule is 
applied to integration. By dividing the z-axis into equally 
spaced elements and evaluating the function at each of these 
points a close approximation of the integral can be obtained 


through the relation 


f(z) da=AZ (£(z,) +4£(z,) +2.£ (zy) +4 (z,) +2£ (25) 


Z, (88) 
+, . 4 WAGE (Zoo) ertaleze)e) 


where an odd number of points must be used [Ref. 12:p. 350]. 
By choosing an endpoint, varying the spacing between points, 
and repeatedly evaluating the integral until no difference is 
observed between successive iterations, an optimal 4z is 
obtained. A close approximation to the semi-infinite integral 
is obtained by using that value of 4z while increasing the 
value of the upper limit on the integral until no difference 


is detected between successive iterations. 
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Figure 3.1: Example of dividing a function into equally 
Spaced segments suitable for use with Simpson’s rule. 
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IV. RESULTS 


A. PREVIOUS RESULTS 

Studies concerning the effects of non-uniform windowing on 
a DFT based receiver in the presence of a Doppler shifted 
Signal have been performed and documented in [1] and [2]. Each 
study draws conclusions concerning the effectiveness of the 
DFT based receiver under different circumstances. The results 
of these studies are included in order to better understand 
the conclusions of this thesis. 

Reference [1] studies the effects of non-uniform windowing 
on a DFT based receiver in a nonfading channel with a Doppler 
shifted signal, where the probability of bit error as a 
function of the signal-to-noise ratio is obtained. It is 
demonstrated that uSing a non-uniform window (specifically, a 
Hamming window) results in a degradation of performance as 
compared with the performance of a uniform window for the 
cases of little of no Doppler shift on the incoming signal. 
However, performance enhancement is obtained as Doppler shift 
increases. For the case of M=8, a Doppler shift of e=.33 is 
sufficient to warrant the use of non-uniform windowing. 
Furthermore, the performance improvement increases as the 
amount of Doppler shift in the signal increases. [Ref. 1: pp. 


2017-2018] 
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The effects of a Rician fading channel on the DFT based 
receiver are investigated in [(2]. Another parameter for 
analysis is the direct-to-diffuse ratio discussed in the 
previous chapter. In a manner similar to that developed in 
(1], data is generated relating the probability of bit error 
to the signal-to-noise ratio for different direct-to-diffuse 
ratios. For large direct-to-diffuse ratios (i.e., conditions 
corresponding to those studied in [1]), the results correlate 
well. However, as the direct-to-diffuse ratio is decreased, 
the non-uniform window gives better performance than a 
rectangular window only when the Doppler shift decreases. The 
presence of a strong diffuse signal, as illustrated by a small 
direct-to-diffuse ratio, limits the benefits of using a 


rectangular window. 


B. EFFECTS OF FREQUENCY-HOPPING AND PARTIAL-BAND INTERFERENCE 

Graphs of the worst case probability of bit error as a 
function of signal-to-jamming noise ratio for fixed signal-to- 
thermal noise ratios are constructed for one, two, three, and 
four hops per symbol. These results are obtained by 
numerically evaluating (40), (48), and (49) for values of y 
ranging from 1.0 to 0.001 and retaining the worst case 
performance for each value of E,,/N; - These results are plotted 
in Figs. 4.1 through 4.4. The signal-to-noise ratio, direct to 
diffuse ratio, and the modulation order are chosen to enable 


direct comparisons with the results obtained in [4] which 
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Figure 4.1: Worst case probability of bit error versus 


Signal-to-jamming noise ratio for L=1, M=4, direct-to- 
diffuse ratio=10, and SNR=13.35 QB. 
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Figure 4.3: Worst case probability of bit error versus 


Signal-to-jamming noise ratio for L=3, M=4, direct-to- 
diffuse ratio=10, and SNR=13.35 dB. 
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Figure 4.4: Worst case probability of bit error versus 
Signal-to-jamming noise ratio for L=4, M=4, direct-to- 
diffuse ratio=10, and SNR=13.35 dB. 
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analyzes a conventional (non-DFT) noise-normalized MFSK 
receiver with quadratic detection. The results obtained here 
for the case of no Doppler shift (e€=0) and a rectangular 
window should be similar to those obtained in [4]. 

Differences in the results obtained in [4] and those 
obtained in this study are primarily due to the facts that the 
equivalent system bandwidths used in [4] and this thesis are 
different and the results in [4] are exact while the results 
in this study represent a union bound. The equivalent system 
bandwidth used in [4] is chosen to be minimal while the 
equivalent system bandwidth used in this work is chosen to 
enable increased frequency spacing between the DFT output bins 
as well as to allow comparisons to be made with the results in 
{1] and [2]. From this, it 1s expected that the data generated 
for this thesis will predict slightly worse performance than 
that predicted an"[(4]. This )aemebserved ianvall casee. 

Just as reported in [4], the performance of the DFT based 
communications system is enhanced as the number of frequency- 
hops per symbol is increased. This is especially true for the 
cases of little or no Doppler shift. For cases of large 
Doppler shift (i.e., e€=0.5), the performance of the receiver 
is slightly poorer as the hop rate is increased when non- 
uniform windowing is used. The error introduced by self- 
induced crosstalk overcomes the benefits usually associated 
with fast frequency-hopping. This is clearly seen in Fig. 4.1 


and Fig. 4.4 for the case of a rectangular window and é«=0.5. 
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For E,/N;=41.0 dB, Pb=4.41x10° for a hop rate of one hop per 
symbol, and Pb=6.65x10° for a hop rate of four hops per 
symbol. The use of a Hamming window in this situation improves 
receiver performance as the number of hops per symbol 
increases, but the benefits are small compared to the 
complexities needed to implement a fast frequency-hopped 
receiver. 

For the case of slow frequency-hopping, the results in 
Fig. 4.1 approach those reported in [2] as E,/N; becomes very 
large. This is expected because the decision statistic 
consists only of one component due to the fact that the 
Carrier frequency is changed once per symbol. Furthermore, 
receiver performance degrades as either the jamming power or 
the Doppler shift increases. As reported in [1] and [2], the 
performance of the Hamming window is worse than that for a 
rectangular window for small Doppler shifts, but outperforms 
the rectangular window for large Doppler shifts. The value of 
€ at which the use of a Hamming window gives’ better 
performance than the use of a rectangular window (the 
crossover fractional frequency offset) is a function of the 
amount of jamming power at the receiver. At E,/N;=1.0 hye Lor 
example, with a rectangular window, Pb=0.413 with e=0 and 
Pb=0.687 with e=0.5; while, with a Hamming window, Pb=0.524 
with e€=0 and Pb=0.668 with e«=0.5. However, at E,/N;=41.0 Ci) 
with a rectangular window, Pb=4.59x10 with «=0 and 


Pb=4.41x10° with ¢«=0.5; while, with a Hamming window, 
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Pb=1¥926x107" wiithe=-OWand Pb=7.87x10~" witWc-085" ror lane 
Doppler shifts, the use of a Hamming window clearly improves 
receiver performance over a broad range of E,/N;- 

This same phenomena exists for the cases of two, three, 
and four hops per symbol. In these cases, however, the larger 
hop rates have little effect on the value of the crossover 
fractional frequency offset. This is illustrated in Figs. 4.5 
through 4.8 which plot the worst case probability of bit error 
as a function of the fractional frequency offset for 
E,/N;=16.0 aB and L=1,2,3 and 4. These results are summarized 
in Fig. 4.9, which plots the crossover fractional frequency 
offset as a function of the hop rate for various values of 
E/N; . As can be seen, the values of the crossover fractional 
frequency offset remain in a small band for each value of 
E,/N;, and increase as the amount of jamming power at the 
receiver increases. 

From Fig. 4.9, it can be seen that the crossover 
fractional frequency offset does not follow a smooth curve as 
the hop rate is varied. In viewing these curves, it must be 
remembered that the results in Figs. 4.1 through 4.4 represent 
the worst case performances of union bound solutions. These 
worst case results are obtained by varying the value of the 
fraction of spectrum jammed (y) and retaining only the worst 
case performance. The value of y giving worst case performance 
at one hop rate is not generally the same value at a different 


hop rate. The fact that only a representative number of 


4 4 


discrete values of y are used to approximate the continuous 


range of possible y's contributes to this situation. 
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Figure 4.5: Worst case probability of bit error versus 
fractional frequency offset for the case of L=1, M=4, 
SNR=13.35 dB, and E,/N,=16.0 dB. 
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Figure 4.6: Worst case probability of bit error versus 
fractional frequency offset for the case of L=2, M=4, 
SNR=13.35 dB, and E,/N,=16.0 GB. 
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Figure 4.7: Worst case probability of bit error versus 
fractional frequency offset for the case of L=3, M=4, 
SNR=13.35 dB, and E,/N;=16.0 dB. 
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Figure 4.8: Worst case probability of bit error versus 
fractional frequency offset for the case of L=4, M=4, 
SNR=13.35 a@B, and E,/N,;=16.0 Ger 
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V. CONCLUSIONS AND RECOMMENDATIONS 


As the frequency-hop rate of the DFT based communications 
system with noise-normalization increases, receiver 
performance improves in the absence of large Doppler shifts. 
However, when large Doppler shifts exist, the errors 
introduced as a result of self-induced crosstalk counteract 
any benefits gained from the use of fast frequency-hopping. 
The advantages gained by using a non-uniform window in the 
presence of large Doppler shifts improve the situation, but 
not enough to warrant the added hardware complexities inherent 
to the application of fast frequency-hopped spread spectrum. 

The amount of Doppler shift necessary to warrant the use 
of a non-uniform window is primarily a function of the jamming 
noise power spectral density. Holding all else constant, as N, 
increases, the amount of Doppler shift necessary to justify 
the use of a non-uniform window increases. At very small 
values of E,/N, the use of a non-uniform window results in 
degraded performance for all but the most extreme Doppler 
shifts. The performance degradation resulting from using a 
non-uniform window in a situation with little Doppler shift, 
however, is small compared to the performance enhancement 
gained by using a non-uniform window when large Doppler shifts 


exist. In general, in the absence of all information regarding 


ol 


the nature of the received signal, the use of a non-uniform 
window is preferred in all situations. 

For all but the most severe conditions of Doppler shift 
and jamming power, the receiver studied in this thesis 
provides reliable communications through the use of fast 
frequency-hopped spread spectrum. Communications reliability 
does not come without a price. The complexities involved in 
implementing a fast frequency-hopped spread spectrum system 
are enormous and are not worth the benefits in all situations. 
As an example, in order to keep the same number of samples per 
frequency-hop for the DFT, the total symbol duration time 
becomes a linear function of the hop rate. In this thesis the 
hop time is determined to be equal to one second. For the case 
of one frequency-hop per symbol, this translates to a symbol 
duration time of one second, but for the case of L frequency 
hops per symbol, this translates to a symbol duration time of 
L seconds. Thus, a faster frequency-hop rate translates to 
slower, but in many situations, more reliable communications. 

Despite its drawbacks, the use of the DFT based receiver 
may be superior in Situations where the content of the message 
is more valuable than the speed of delivery. By employing the 
DFT based receiver with fast-frequency hopping, a system can 
be developed which has the potential to deliver reliable 
communications with considerable anti-jam protection. 
Situations relying on this type of communications protection 


are common in the military environment. 
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In order to gain deeper insight into the relationships 
between the fractional frequency offset, the modulation order, 
the hop rate, and the signal energy-to-jamming noise power 
spectral density ratio, further study should be conducted. 
This study should emphasize the changes in the crossover 
fractional frequency offset observed while varying the key 
parameters listed above. Also, the calculation of the exact 
solution is possible, and should be utilized in any further 
study. Due to limited computational resources, this was not 
possible in this thesis. 

Finally, the work presented here is theoretical in nature 
and should be tested with actual hardware. The assumptions 
made in this analysis were selected to ease the burden of 
computational effort needed to solve the problem. They may not 
correspond to - actual conditions needed for optimal 
communications in a real world system. Only building this 
receiver and testing it under realistic conditions will prove 
or disprove the usefulness of the DFT based, noise-normalized 


receiver. 
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